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ABSTRACT 

We study the excitation and damping of transverse oscillations in a multi- 
stranded model of a straight line-tied coronal loop. The transverse geometry of 
our equilibrium configuration is quite irregular and more realistic than the usual 
cylindrical loop model. By numerically solving the time-dependent ideal mag- 
netohydrodynamic equations in two dimensions we show how the global motion 
of the whole bundle of strands, excited by an external disturbance, is converted 
into localized Alfvenic motions due to the process of resonant absorption. This 
process produces the attenuation of the transverse oscillations. At any location 
in the structure two dominant frequencies are found, the frequency of the global 
mode, or quasi-mode, and the local Alfven frequency. We find that the mech- 
anism of mode conversion, due to the coupling between fast and Alfven waves, 
is not compromised by the complicated geometry of the model. We also show 
that it is possible to have energy conversion not only at the external edge of the 
composite loop but also inside the structure. The implications of these results 
and their relationship with the observations are discussed. 

Subject headings: MHD — Sun: corona — Sun: magnetic fields — waves 

1. Introduction 

The resonant coupling of fast magnetohydrodynamic (MHD) waves and Alfven waves 
has been an active research topic for many years. The most important applications are to 
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laboratory plasmas, the solar corona, and pulsations in the Earth's magnetosphere. In the 



didate to explain the damping of transverse coronal loop oscillations (see 


HoUwee & Yang 


1988; 


Goossens et al. 


1992; 


Ruderman & Roberts 


2002; 


Terradas et al. 


2006a 


). The attenu- 



ation of the loop oscillations, interpreted as the kink modes, is due to the energy conversion 
from the global large scale motions to localized Alfven modes. 

The time-dependent resonant absorption mechanism has been traditionally analyzed 
for both driven and initial value problems using different equilibrium mo dels. For the 
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Ofman &: Davila 



in cy l indrical tubes (IGrossmann fc Smith 



1995 



driven problem this mechanism ha s been studied in smooth transition layers (iHoUweg 
Davilalll987l: iHoUweg fc Yand|l988[). in slab naode. 



lOfman et a. 



1991; Sakurai et al. 1991; Goossens et al 
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mode 
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(Steinolfson & Davila 


1993; 


Ofman et al. 


1995 




De Groof & Goossens 


200C 




2002). and 



Poedts et al 



1989 



199G 



1987 



Poedts fc Kerner 



For the in i tial value problem sim ilar stud- 



ies have also focused on smooth interfaces (lonson 



1978 



Lee fc RobertsI 1986), on slabs 



he ef f ects of backg r ound fl ows have been i nvesti g ated bvlHoUweg et al 



( Steinolfson fc Davilalll993 ). and cylinders ( Ruderman fc Robertslboo"? Terra.das et al.ll2006a ). 



(Il996l ): iTirrv et all Jl998h ; lAndries et~aD J200oh : kndries fc GoossensI ( l200lh . 



1990l ): lErdeyi fc Goossens 



The tin ae-dependent results hav e been complemented by eigenmode calculations (see 
for example iTirry fc Goossenslll996l ). The coupling of the fast kink mode with the Alfven 
continuum produces the damping of the oscillations. In this case the global mode is usually 
called the quasi-mode and has a complex frequency. The theoretical models that have 
been used to calculate the quasi-modes are so far simple but necessary to understand the 
main properties of resonant absorptio n. Recently, new effects have be en in vestigated. The 



curvatu re effect has been studied by IVan Doorsselaere et al.l (l2004bl ) and 
( 2006br) . and t he influence of s tratifi cation along the loop has been analyzed by 



erradas et al. 



Andries et al. 



kood ) and by lArregui et all J2005h. The effect of non-circularity of the tube cross-section 



has been analyzed by iRudermaru (12003), while the effect of the internal structure in loops 



using a slab model, has been studied by lArregui et al.l (120071 ). The overall conclusion of these 
recent investigations is that these new ingredients do not significantly change the damping 
per period of the oscillations. This suggests that resonant absorption is a robust mechanism 
and that its efficiency is not easily affected by the considered effects. 

Nevertheless, even the most recent models are still too simple compared with the real 
conditions in coronal loops. For example, there are no physical reasons to think that the 
cross-section of coronal loops is perfectly cylindrical. Moreover, there is observational infor- 
mation suggesting that loops are not monolithic (as they are usually modeled) but that they 
are formed by bundles of individual strands considered as mini-loops for which the heating 
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/uuol ) . i ni s view is supported by some autnors (iMartens e 
by others ( Aschwanderj 2005 : Aschwanden &: Nightingalel 



and p lasma properties are approximately unifo rm in the transverse direction (Schnaelz et al. 
20051 ) . Thi s view is supported by some authors (IMartens et al.ll2002l : lKlimchukll2006l ) but not 



20051 ). Thus, one of the questions 



that arise from these observations is how the internal structure of coronal loops can modify 
the mechanism of resonant absorption. 

The purpose of this paper is to study the mechanism of resonant absorption and its 
efficiency in a rather arbitrary two-dimensional distribution of plasma, and to analyze the 
implications on the damping of transversal coronal loop oscillations in a multi-stranded loop 
model. Instead of calculating the eigenmodes of the structure, which is difficult due to the 
complicated geometry, we investigate the dynamical response of the bundle of individual 
strands to an initial perturbation. Thus, we study the initial value problem of the excitation 
of the composite loop by solving the time-dependent two-dimensional ideal MHD equations. 
The equations are solved numerically using an appropriate code. From the simulations we 
study the motions in the loop structure and analyze how the energy from the global mode 
is converted into localized Alfvenic motions in the inhomogeneous regions of the loop. 

This paper is organized as follows. In § [2] the multi-stranded loop model is described 
and the basic MHD equations are given. The numerical method used to solve the time- 
dependent equations is explained in § [31 In § H] the results of the time-dependent problem 
for a single cylindrical inhomogeneous loop are compared with the eigenmode calculations. 
In § O the time-dependent problem for the multi-stranded loop is considered and several 
features of the resonant absorption mechanism are discussed, in particular the behavior in 
the inhomogeneous layer, and the energetics of the problem is investigated. Finally, in § [6] 
the main conclusions are drawn. 



2. Loop Model and Governing MHD Equations 

The equilibrium magnetic field is straight, uniform, and pointing in the 2;— direction, 
B = i?o For applications to the solar corona, it is a good approximation to consider that 
the magnetic pressure dominates the gas pressure. This zero— /3 approximation allows us to 
choose an arbitrary density profile. As a model of a bundle of loops we use a superposition 
of tiny, parallel cylinders with different radii and densities. In Cartesian coordinates, the 
cross-section of the density of each individual strand is assumed to have the following form. 
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where poi is the maximum density of the strand, Xi and yi the position of the strand axis, 
and tti the strand radius. The density of the multi-stranded model is defined as 



N 



Po{x, y) = ^ Piix, y) + Pex, (2) 



1=1 



Pex being the density of the external medium (assumed constant). Note that the density of 
a single uniform cylindrical loop, pin, with radius R that has the same mass as the multi- 
stranded loop is simply, 

1 ^ 

= X] + P^""- 

i=l 

In Figure [T]the two-dimensional distribution of the density (the cross section of the composite 
loop) is plotted for a particular configuration based on equations ([1]) and ([2]). For this 
particular configuration the loop is composed of ten individual strands with their axis located 
at the following Xi and yi coordinates: [0, 0, -0.65, 0.65, 0.5, -0.5, -0.5, 0.5, 0.2, -0.2] and 
[-0.45, 0.65, 0, 0, 0.5, 0.5, -0.5, -0.5, 0, -0.2]. The radii of the strands (ai/R) are [0.2, 0.275, 
0.25, 0.3, 0.2, 0.25, 0.2, 0.25, 0.25, 0.25], and the densities (poi/poo) are [0.8, 0.6, 0.4, 0.3, 0.5, 
0.6, 0.4, 0.6, 0.4, 0.3]. The external density is pex = 1/3 poo. Poo being a reference density 
inside the loop. 

We see that the loop density has an inhomogeneous distribution with quite an irregular 
cross section and irregular boundary. This model has a complex geometry compared with 
the usual cylindrical or the el hptical tubes usually used to study loop o scillations. In this so- 



called "spaghetti model" (see lBogdan fc Foxlll99ll : lKeppens et al.lll994j ) the distance between 



the strands is quite small and a strong dynamical interaction between them is expected. This 
kind of model has been previously considered in the context of scattering and absorption of 
sound waves in composite sunspot models. 

To study small amplitude perturbations in this equilibrium we use the linearized ideal 
MHD equations. The equilibrium depends on x and y but it is independent of the longitudinal 
2;— coordinate. For this reason, we Fourier analyze in this direction assuming a dependence 
of the form e~'^^^^, being the vertical wavenumber. We concentrate on the fundamental 
mode, with kz = tt/L, L being the length of the loop (here we use L = 20R). Under these 
assumptions the MHD equations are, 

Po^ H ^kzbx + = , (4) 



ai p \ ox 
dvy Bq [ dhz , ^ 



- 5 - 




-1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 
x/R 



Fig. 1. — Cross section of the density of the multi-stranded loop model. The lengths are normalized 
to i?, the reference radius of an individual homogeneous loop. The contour lines represent curves of 
constant Alfven frequency, given by equation (jlOp (they also represent contours of constant density) . 
The thick line corresponds to the Alfven frequency that matches the frequency of the global mode. 
The labels 1, 2, and 3 indicate the locations analyzed in Figures El [U andEJ respectively. 



9hx ., 

— + Boik^v^ 
at 

— + Boik^Vy 
dh. ^ f dvr. du 



dt 



dx ~'~ dy 








(6) 
(7) 



where v = {vx,Vy,0) is the velocity and b = {bx,by,bz) is the perturbed magnetic field. 
These equations can be rewritten as a system of five equations with five real variables using 
the transformations ^ ib^ and by iby (these variables are shifted 7r/2 in the z-direction 
with respect to v^, Vy, and b^)- The equations are complemented with an initial perturbation 
located in the external medium. For simplicity we assume that it has the following form 



Vy{x,y) = Vyo exp 



(9) 



where Vyo is the amplitude of the perturbation (this parameter is arbitrary since we are in the 
linear regime), yo is the position of the center of the disturbance, and w is its width (here we 
use yo = 3R and w = R, which means that the perturbation is located at a distance 3-R from 
the center of the loop). All other MHD variables are initially set to zero. This perturbation is 
a planar pulse which produces the excitation of fast MHD waves that propagate and interact 
with the loop structure. 
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In the excitation of fast MHD waves and specially Alfven waves an important quantity 
in our model is the local Alfven frequency, 

ujAix,y) = k^vx{x,y) = ° =, (10) 

va being the Alfven speed. The Alfven frequency distribution in our equilibrium depends 
on X and y in a quite complicated manner, as we can see in Figure [H (see the contours of 
constant uja represented with continuous lines). The fact that uja varies with position plays 
a key role for the resonant conversion of wave energy from global large scale motions to 
localized Alfven modes. 



3. Method 

From the numerical point of view the problem that we are studying has basically two 
difficulties. The first one is related with the small spatial scales that are generated in 
the inhomogeneous layers, so we must make sure that we are properly resolving these lay- 
ers. This point is crucial to have the correct energy conversion rates and in consequence 
the right damping times of transversal oscillations. Due to the phase mixing process (see 



Heyvaerts &: Priestlll983l ). specially important at the inhomogeneous layers, the maximum 



time for which the simulations can be run is determined by the phase mixing le ngth, defined 



as the length over which the phase of neighboring Alfven waves differs by 27r ( iMann et al. 



19951 : IWright fc RickardlllQQSf ). 



(u) 

Therefore, the typical spatial lengths decrease quickly with time. The shortest wavelength 
that can be resolved numerically with a uniform grid spacing of Ax is A = 2 Ax. We can run 
the simulations until the phase mixing length is equal to the shortest wavelength (Lph = A). 
Thus, the maximum simulation time is simply 

tmax = vr/ (^Ax . (12) 

We have checked that this condition is not violated in the simulations (note that it also has 
to be satisfied in the y— coordinate). 

The second difficulty is the effect of the boundaries on the loop dynamics. Even applying 
transparent boundary conditions (zero order extrapolation) at the limits of the computational 
domain, it is important to locate the boundaries far enough from the loop. Note that a 
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trapped mode in the loop has always an evanescent part in the external medium which 
might be affected by the conditions that we apply at the boundaries. 

To solve these difficulties we have used a resolution sufficiently high to resolve the 
different scales (a grid of 4000 x 4000 points is sufficient) and to avoid significant reflections 
we have located the domain limits far enough from the loop, at xb = ±8-R, i/b = ±8-R. 
Nevertheless, in order to better interpret and visualize the results the plots are displayed in 
a smaller spatial domain ([—1.6/2, 1.6-R] x [— 1.6-R, 1.6-R]). 

To numerically solve the time-dependent MHD equations, equati ons dH) — (El), to gether 



I20 



with the initial condition, equation ([9]), we use the code CLAWPACK (lLevequell2004l ). This 
code implements a wide class of high-resolution finite volume methods for solving linear or 
nonlinear hyperbolic problems. Due to resolution requirements we have run the parallelized 
version of the code in a cluster of computers. 



4. Test case: resonant absorption in an inhomogeneous cylindrical loop 



Due to the different scales involved in the problem the modeling of resonant absorp- 
tion is challenging from the numerical point of view. For this reason, before the "spaghetti 
model" is studied, it is necessary to check that the numerical method we are using is appro- 
priate to study this problem. As a test we consider a single inhomogeneous cylindrical loop 
in two dimensions. To facilitate comparison of the time-de pendent results with previou s 
eigenmode studies we choose t h e same density profile as in iRuderman fc Roberts! (120021 ): 
Van Doorsselaere et al.l (l2004al ): iTerradas et al.l (l2006al ). i.e., a sinusoidal variation in den- 
sity across the non-uniform layer. For such a configuration and for the m = 1 mode (the 
kink mode), the damping per period (the damping time, Td, over the period, P) in the limit 
of thin tube and thin boundary is 

Td 2 i? Pin + p 



(13) 



TT / Pin - Pe 

where I is the thickness of the non-uniform layer and R is t he loop radius. For thick layers 
the eigenvalue proble m has to be solved numerically (see IVan Doorsselaere et al.l l2004al : 
Terradas et al.ll2006al ). In our composed loop model there is a wide range of thicknesses, 
and for this reason we ha ve calculated the eigenmodes of the cylindrical loop numerically 



see 



Terradas et al.ll2006al . for details about the method used to perform these calculations). 



In Figure [21 (top panel) the damping per period for a single cylindrical loop calculated from 
equation (1T31) and the eigenmode calculations are represented as a function of the thickness 
of the layer. The small differences between the two curves are simply due to the fact that 
equation (1T31) is inaccurate for thick layers (see also I Van Doorsselaere et al.ll2004al ). 
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Fig. 2. — Top panel: Damping per period as a function of the thickness of the non-uniform layer 
for a single cylindrical loop (pex = 1/3 /Om, kzR = 7r/20). The dashed line is the damping per 
period calculated using equat ion (fTSll while the con tinuous line is calculated with the eigenvalue 



problem for thick layers (see iTerradas et alj l2006al ) . The circles represent the values calculated 
from the time-dependent two-dimensional simulations. Bottom panel: Plot of Vy at the center of 
the cylindrical loop as a function of time for the case l/R = 0.6. This result has been obtained by 
solving the time-dependent equations. The dashed line is a fit of the form exp(— t/r^). The time is 
normalized to the Alfven transit time, ta = VAo/R, where vao = -Bo/VWin- 

To compare the results of these eigenmode calculations with the time-dependent problem 
we need to place an initial perturbation in the system that basically excites the kink mode 
(the eigenmode calculations that we want to compare with are for this mode, m = 1). 
The planar pulse perturbation given by equation (Q is the kind of disturbance with such 
a property because an external perturbation hardly excites high order har monics such as 



the flu ting modes, and the sausage modes are leaky in our conflguration (see iTerradas et al. 



2007bl ). We have run the code with this perturbation and have studied the subsequent 



evolution. Figure [2] (bottom panel) shows the plot of Vy at the c enter of the loop as a 



function of time for a particular thickness of the layer (see Fig. 4 in ITerradas et al.ll2006al . 
for the analogous simulation in ID). After a short transient the loop oscillates with the quasi- 
mode period and the amplitude of the mode is attenuated due to the energy conversion in 
the inhomogeneous layer. From the results of the simulations we calculate the period, P, 
by performing a periodogram. The damping time, t^, is numerically estimated by fltting an 
exponential function of the form exp(— t/r^) to the envelope of the curve. We have performed 
several simulations varying the width of the layer and the numerically obtained damping per 
periods are represented with circles in Figure [2] (top panel) as a function of the thickness 
of the layer. We flnd the expected dependence with l/R and we see that the deviations 
from the eigenmode calculations are quite small. It is important to remark that, although 
we model a cylindrical tube, in these simulations we use Cartesian coordinates and we still 
obtain the right damping rates. 
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In conclusion, with this simple numerical experiment we demonstrate that we are using 
a reliable tool to study resonant absorption. Note that contrary to some precedent works, 
where a large amount of resistivity is used mainly to avoid numerical problems, here we are 
in ideal MHD conditions. Hence, the system dynamics is not dominated by the resistive 
regime, in this case the behavior in Figure [2] (top panel) would be completely different, the 
damping time being independent oil/R (see for example Fig. 2 in lTerradas et al.ll2006al ). 



5. Results: multi-stranded model 

We now use the multi-stranded model shown in Figure [T] and study how the system 
evolves due to a planar pulse perturbation (eq. [9]). The initial pulse produces a displacement 
of the whole ensemble of strands basically in the y— direction, i.e. the direction in which 
the initial perturbation propagates. The initial stage of the evolution of the bundle of 
loops is dominated by a complicated set of internal reflections of the wavefront between the 
different strands. During this transient several wavefronts propagating from the bundles 
into the external medium are found. They correspond to the emission of the leaky modes 



representirig fast radiating MHD waves with short periods and fast attenuation (see ICally 



19861 . l2003l : iTerradas et al.ll2007al . for the analysis of such modes in a single cylindrical loop). 
After the transient the system periodically oscillates with certain frequencies. Hereafter, we 
concentrate on this stage of the dynamics. 



5.1. Frequency analysis 



We first investigate the characteristic frequencies of the system. An analysis of the 
results of the simulations shows that inside the loop there are basically two dominant fre- 
quencies at each point. One is the collective frequency of the bundle of loops (different from 
the kink frequency of the individual strands) and the other is the local Alfven frequency, 
t^A- The collective frequency is the result of the excitation of the global mode of the system 
which represents the emergent behavior of the entire loop. On the other hand, the local 
Alfven frequency is due to the excitation of the Alfven continuum modes. As an example, in 
Figure El the Vy velocity component as a function of time (top panel) and its power spectrum 
(bottom panel) are represented at an interior point (see label 1 in Fig.[T]). The two dominant 
frequencies, the local Alfven frequency and the global frequency (with the largest power) are 
clearly identified in the power spectrum. Note the strong decrease of the amplitude with 
time. This attenuation is found almost everywhere at the interior points and suggests a 
collective motion of the bundle of tubes. We have also represented in Figure [3] (top panel) 
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the other component of the velocity, Vx- Its amphtude is smaller that the Vy component and 
it does not show attenuation with time. 




Fig. 3. — Top panel: Time evolution of the (dashed line) and Vy (continuous line) velocity 
components at the interior point x = 0, y = 0.5R (see label 1 in Fig[T]). Bottom panel: Power 
spectrum of the Vy velocity component. The peak with the largest power corresponds to the 
global mode while the weaker peak is due to the excitation of the local Alfven mode. The local 
Alfven frequency, calculated from the function Lii\(x,y), at the coordinates of the previous point is 
represented with dot-dashed lines and agrees with the location of the weak peak. 




Fig. 4. — Top panel: Time evolution of the Vx (dashed line) and Vy (continuous line) component 
at a position near the loop edge, x = 0.9R, y = (see label 2 in Fig{T]). Bottom panel: Power 
spectrum of the Vy velocity component. The local Alfven frequency is represented with dot-dashed 
lines. 



There are particular locations in the structure where the behavior of the signal is com- 
pletely different. An example is shown in Figure H] which corresponds to a position situated 
near the loop edge (see label 2 in Fig. [1]). Now, instead of a damped oscihation, the ampli- 
tude of the Vy velocity component grows with time while the Vx component decreases. There 
is just a single peak in the power spectrum because the local and the global frequencies coin- 
cide. As we will show in the following sections this is the place where the energy conversion 
takes place. 
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Fig. 5. — Top panel: Time evolution of the Vx (dashed line) and Vy (continuous line) component 
at a position near the loop edge, x = lAR, y = (see label 3 in Fig{T]). Bottom panel: Power 
spectrum of the Vy velocity component. The local Alfven frequency, i.e. for this particular point 
the external Alfven frequency, is represented with dot-dashed lines. 



Outside the structure, oscillations tending to the Alfven frequency are found everywhere, 
see for example the velocity represented in Figure [5] (see label 3 in Fig. [1]). These oscillations 
are basically produced by the wake of the initial fast MHD perturbation and have a large 
amplitude compared with the amplitude inside the loop (see Fig. [3D. This wake os cillates 



(in the limit of large times) at the local Alfven frequency (see iTerradas et al.l l2005l ) which 
is precisely where the dominant peak is found in the power spectrum. However, it is worth 
noticing that near the loop edge, but in the external medium, the quasi-mode frequency 
is also present. This is not surprising, since the quasi-mode has an evanescent tail in the 
external medium. 

We have also performed the same analysis but for the bz component, which is propor- 
tional to the total pressure perturbation (Pt = bzBo/fi). We have found large power at 
the global mode frequency and very weak power around the local Alfven frequency. This is 
because pure Alfven modes do not perturb the total pressure and, although in our config- 
uration Alfven modes are coupled with fast modes, they still keep a strong incompressible 
character. 



5.2. Velocity field 

To explain the different behavior of the velocity at different locations in the structure 
we need to understand the evolution of the whole system instead of looking at individual 
points. In Figure E] the velocity components v^, Vy are plotted at t = 40 r^. Inside the loop 
the Vx component is quite uniformly distributed while the Vy component has a complicated 
spatial structure which basically coincides with the spatial distribution of the strands (see 
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the contours). At this particular instant the higher the density of the strand the smaller the 
amplitude of the Vy component. On the other hand, the velocity components outside the 
loop have quite a smooth spatial distribution. The motion of the whole structure is clearer 
if we represent the velocity field constructed using the and Vy components. The result is 
shown in the top panel of Figure [71 We clearly see that all the strands are moving in the 
negative ?/— direction and that the external medium reacts in a very organized way to the 
displacement of the bundle of loops. The plasma at the bottom of the loop {y ~ —R) is 
pushed sideways while the plasma at the top of the structure {y ~ 0.8 i?) tends to fill the 
region that has displaced downwards. This motion is the equivalent of the kink mode in a 
cylindrical tube. 




-a.E-02 -4.E-02 O.E + 00 4.E-02 B.E-02 



Fig. 6. — Spatial distribution of the Vx velocity component {top panel), Vy velocity component 
{bottom panel) at t = 40 r^. The corresponding velocity field is plotted in the top panel of Figure [71 
Contours of constant Alfven frequency are represented with continuous lines. 

At later times (see the middle and the bottom panels in Fig. [7]) large amplitude velocities 
develop specially near the external edge of the composite loop. This is the consequence of 
the energy conversion between the global mode and the Alfven modes. Due to this process 
the global oscillation gradually losses its energy and its amplitude is attenuated in time (see 
for example the length of the arrows at the centers of the individual strands in Fig. [71 at 
t = 160 r4). This explains the attenuation of Vy found in Figure [31 On the other hand, the 
amplitude of the Alfven modes increases at the resonant layers where the energy conversion 
takes place (see the large arrows at the loop boundaries in Fig. [71 bottom panel). We already 
found this behavior in Figure [H 



- 13 - 



V^A= 40 



1.5 - / 

- , / 

: 


////y^.^^. . 


J, 


i.5 - 

1.0- 


o.s- 






0.5- 


X 0.0- j 






X 0.0- 


-0.0- 






-0.5- 










-r.o- 






-1.0- 


-1.5- 






-1.5- 



-1.5 -1.0 



-0.5 0.0 0.5 1.0 1.5 
x/R 



-1.5 -1.0 -0.5 



x/R 



llllll 

0.5 1.0 1.5 



t/T,= 160 



1 : 



m\ t 



-1.5 -1.0 -0.5 0.0 5 1.0 1.5 
x/R 



Fig. 7. — Time evolution of the velocity field. The loop initially oscillates in the y direction 
with the global mode. Due to the mode conversion motions localized on the magnetic surfaces 
develop (see the contours of constant Alfven frequency represented with continuous lines). The 
region marked with a square box in the bottom panel is plotted in Figure [8l This figure is also 
available as an mpeg animation in the electronic edition of the Astrophysical Journal. 



Since the Alfven frequency changes with position the Alfvenic motions get out of phase 
very quickly due to phase mixing. The consequence of this process is visible at the external 
boundaries where strong shear motions develop. To see these motions in detail we have 
concentrated on the small region marked with a square box in Figure [7] bottom panel. The 
result is plotted in Figure [HI We see that the field is organized and that near the edge of 
the loop the velocity vectors are aligned with the contours of constant Alfven frequency (see 
the continuous lines). It is also clear that on neighboring magnetic surfaces motions are 
basically in opposite directions. This is even clearer in Figure [9] where we have plotted a 
cut of Vy at y = (see the continuous line). This figure also shows how this component of 
the velocity field (which is mainly parallel to the magnetic surfaces around x = R) evolves 
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Fig. 8. — Detail of tlie velocity field at t = 160 r^, see square box in the bottom panel of Figure[71 



with time. As expected from the phase mixing process, the typical spatial scales of the 
shear motions decrease with time. The amplitude of the velocity decreases in the internal 
part of the structure while it develops short wavelengths around the resonant positions. 
We find that the typical wavelengths agree with the phase mixing lengths, Lph, predicted by 
equation (ITTll (compare the lengths of the lines in the legend of the plot with the wavelengths 
of Vy aX X = R). 




Fig. 9. — Plot of Vy as a function of a; (y = 0) at t = 40 (dot-dashed line), t = 100 ta (dashed 
line), and t = 160 (continuous line). The lengths of the lines plotted in the legend correspond 
to the phase mixing lengths, Lph, predicted by equation (fTTj) at the position x = (y = 0). 



At this point it is worth noticing that due to mode coupling Alfven modes are excited 
everywhere in the structure even at locations far from the resonance. The excitation of 
these modes is already known, specially in driven problems (see for example iMann et al. 



- 15 - 



19951 : iTirry et al.l 119971 : iGoossens fc De Grooil l200ll ). but there are still several questions 



about these modes that need to be addressed, as for example the amount of energy that is 
deposited in the Alfven modes. In any case, far from the resonance, the amplitude of these 
modes is smaller that the amplitude of the global mode (on average). We see in Figure [3] 
(bottom panel) that for this internal position the amplitude of the global oscillation, in the 
Vy component, is around 1.5 times the amplitude of the Alfven modes (the power of the 
peaks, proportional to the square of the amplitude, are around 24 and 10). A more detailed 
analysis, out of the scope of this work, is required. 



5.3. Energy distribution 



The process of mode conversion is even clearer if we study the energetics of the system. 
For this reason we now focus on the wave energy density. Since in our model we adopt the 
zero—/? approximation there are only contributions from the kinetic and magnetic energy to 
the total wave energy. 



E 



1 



(14) 



This quantity is represented in Figure [TO] at different times (same as in Fig. [7]). The plots 
show that the system evolves from a situation where the energy is more or less uniformly 
distributed to a state in which it is concentrated around the resonant layers. The wave 
energy is spatially distributed in such a way that it follows the irregular geometry of the 
loop edge (see the bright areas on the left and right loop boundaries). Out of the resonances 
the energy decreases with time, see the interior points, except for the area around x = 0.2R, 
y = —0.2R where we can find some energy for large times. It should be also noted that the 
system evolves so that the resonance energy width decreases wi th time. Thi s behavior of the 
wave energy at the resonance has been described in detail by iMami et al.l (I19951) (see als o 
the equivalent results in a cylindrical inhomogeneous loop. Fig. 8, in lTerradas et al.ll2006al ). 
To see these effects more clearly we have concentrated on a slice at y = (see Fig. [TTj) and 
have plotted the wave energy density at different times. This figure shows that the energy 
width decreases with time at the external boundaries, basically at x = R and x = —R. 
We also see the peak around x = 0.3R, and the energy decrease with time in the range 
-0.8 < x/R < 0.1 and 0.5 < x/R < 0.7. 
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Fig. 10. — Time evolution of the energy distribution for the same simulation as in Figure [71 The 
three square boxes labeled with a, b, and c in the bottom panel mark the regions where the velocity 
averages have been computed to determine the damping time (see Fig. [T2]l . This figure is also 
available as an mpeg animation in the electronic edition of the Astrophysical Journal. 



5.4. Resonant magnetic surfaces 



It is quite evident from the velocity and energy analysis where the energy conversion 
takes place. However, we can be more precise in the determination of the resonant magnetic 
surface or surfaces. Once we know the frequency of the global mode the location of the 
resonances is basically where this frequency coincides with the local Alfven frequency. In 
Figure [1] contours of the Alfven frequency, uja, are represented together with the global 
frequency, see the contours with thick lines. In this plot we find that there are two resonant 
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Fig. 11. — Wave energy density, calculated using equation (jl4p . as a function of x {y = 0) for three 
different times (same as in Fig. \T0\ the same notation as in Fig. [9] has been used for the lines). The 
energy concentrates around the external edges of the loop (some part also concentrates inside the 
structure, around x = 0.3R). The width of the energy peaks around the resonant layers decreases 
with time. Note also the decrease of energy in the internal part of the loop. 

magnetic surfaces. One surface is located at the external edge of the composite loop and the 
other resonant magnetic surface is located inside the structure (see the small hole around 
X = 0.2R, y = —0.2R). It is precisely at these locations where the energy maps show 
enhanced energy, see for example Figure [TUlbottom panel, and also where the highest velocity 
amplitudes are found, see Figure [71 Thus, for this particular multi-stranded configuration, 
resonant absorption not only takes place at the external boundaries, but some energy is also 
deposited in the internal part of the composite structure. 



Due to the simultaneous excitation of the quasi-mode and the local Alfven modes, the 
determination of the damping rates of the loop oscillation can be quite difficult (see for 
example the time series in Fig. [3]). However, there is a simple way to estimate the damping 
time. We have taken uniformly distributed regions inside the loop and have added up the 
velocity field at a given instant. This basically averages out the contribution of the local 
Alfven modes since these modes develop small spatial scales. At the same time the average 
enhances the coherent behavior of the global mode, having a large spatial length. In Figure [12] 
we have represented the averaged signal as a function of time for three regions inside the 
loop (see the boxes labeled with letters a, b, and c in Figure [10] bottom panel). We see that 
we get almost the same time dependence for the different regions in the loop, indicating that 
the global mode is dominant everywhere inside the loop. Although the structure is quite 
inhomogeneous and irregular the damping time is basically the same everywhere. Now from 



5.5. Damping time 
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the averaged signal we can calculate the damping time using the same method as in the 
comparison with the cylindrical tube, i.e. by fitting an exponentially decaying function. We 
finally find that the damping per period is Td/P = 2.3. 




50 100 150 200 



Fig. 12. — Time evolution of tlie spatially averaged Vy component as a function of time. The 
continuous, dot-dashed, and dotted lines correspond to the (normalized) averages in the square 
regions marked with labels a, 6, and c in Figure [10] (bottom panel), respectively. The profile of the 
three signals is almost the same since the loop is oscillating with the global mode. The damping 
per period is r^jP = 2.3. 



6. Discussion and Conclusions 



We have studied the mechanism of resonant absorption by solving the time-dependent 
problem of the excitation of oscillations in a complicated multi-stranded coronal loop. We 
have shown that the mode conversion takes place in quite irregular geometries like the one 
studied in this paper and that regular magnetic surfaces (considered in previous works) are 
not necessary for this mechanism to work efficiently. This suggests that resonant absorption 
is quite a robust damping mechanism. Although we have analyzed a particular system, the 
behavior found in the present equilibrium is expected to be quite generic of inhomogeneous 
plasma configurations. Since inhomogeneity is certainly present in coronal loops, the resonant 
coupling between fast and Alfven modes can hardly be avoided in a real situation. For this 
reason, resonant absorption seems to be quite a natural damping mechanism. 

The fact that the loop is not monolithic does not affect much the global oscillatory 
properties. Although the loop is composed of different strands we find a dominant frequency 
almost everywhere in the structure, indicating a global motion. It is interesting to com- 
pare the frequency of this global mode or quasi-mode with the frequency of the equivalent 
homogeneous loop with radius R and the same mass. Using equation ([3]) we find that the 
equivalent cylindrical loop should have an internal density = 0.67poo- The kink mode 



- 19 - 



period for a cylindrical tube with this internal density is P = 28.4 ta. This value is in good 
agreement with the period of the quasi-mode estimated from the numerical simulations of the 
multi-stranded model, which is around P = 28.0 ta- This indicates th at the internal struc - 
ture does not change much the global behavior of the loop (see also lArregui et al.l 120071 ). 
Nevertheless, it may have some effect on the location of the energy deposition. We have 
found that in our model there is also energy conversion inside the loop, although it is small 
compared with the energy at resonances in the external layers. 

Although at any time and at any position both the quasi-mode and the local Alfven 
modes are excited, we have been able to estimate the damping time of the quasi-mode by 
performing averages of the dominant velocity component in different regions of the loop. The 
average eliminates the local Alfven modes and enhances the global mode and is basically what 
the observations provide, a sort of integration along the line of sight (but in the displacement 
instead of the velocity). 

On the other hand, it must be noted that a system of N tubes, like the one studied in this 
paper, is expec ted to have a large number of eigenmodes. For example, in a configuration with 
just two loops iLuna et al.l (120071 ) have shown that there are four kink-like eigenmodes, and 
that an external disturbance in general excites these four modes. However, these authors have 
shown that since their frequencies are very similar it is very difficult to distinguish between 
the different eigenmodes in a time-dependent study. For this reason, the interpretation of the 
global mode found in our multi-stranded loop as a sum of different eigenmodes with similar 
frequencies and similar damping times needs to be considered. This could explain the small 
differences in the period and damping time of the velocity averages in different regions of 
the structure. A detailed analysis of this issue is out of the scope of this paper but it is 
clear that the calculation of the eigenmodes of complicated configurations is very important. 
In this regard, analytical studies b ased on scattering theory will allow us to calculate the 
eigenmodes of a system of tubes ( jLuna et al.ll2008l . in preparation) and to make progress 
in this direction. 



Finally, it must be noted that we have concentrated on the linear regime. If the ampli- 
tude of the oscillations becomes large in the inhomogeneous layer due to mode conversion, 
non-linear terms might be important and the efficiency of resonant absorption can be altered. 
However, diffusion or viscosity may prevent the development of large amplitudes around the 
resonant layer. In addition, there are some results that indicate that in the non-linear regime 
the heating at the resonant layers may produce s ignificant changes in the equilibrium config- 



uration (see lOfman et al.lll998l : iKlimchukI |2006| ) . Hence, the initial value problem needs to 
be studied using the full non-linear equations, which basically means that, since the Fourier 
analysis is no longer possible in the direction, the problem is three-dimensional. Due 
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to the grid resolution requirements in the layers the three-dimensional problem has a high 
computational cost. The numerical study presented here for the two-dimensional problem 
is a preliminary step to investigate resonant absorption in more realistic three-dimensional 
models including, for example, twisted or tangled magnetic fields. 
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